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Abstract 

We study the properties of low- lying Dirac modes in quenched compact QED at (3 = 1.01, 
employing 12 3 x Nt (Nt = 4, 6, 8, 10, 12) lattices and the overlap formalism for the fermion 
action. We pay attention to the spatial distributions of low-lying Dirac modes below and 
above the "phase transition temperature" T c . Near-zero modes are found to have universal 
anti-correlations with monopole currents, and are found to lose their temporal structures 
above T c exhibiting stronger spatial localization properties. We also study the nearest- 
neighbor level spacing distribution of Dirac eigenvalues and find a Wigner-Poisson transition. 

1 Introduction 

In the previous paper [T] , the spatial correlations between monopoles [5] and low- lying Dirac 
eigenmodes (eigenfunctions of the Dirac operator) were studied in 4D quenched compact 
QED, and it was found that there exist universal anti-correlations between them below and 
above the critical coupling (3 C . The clear anti-correlation between Dirac eigenfunctions and 
monopoles implies that the level dynamics of Dirac eigenvalues, which is responsible for the 
chiral phase transition via Banks-Casher relation [3J, is controlled by monopole configurations 
in 4D Euclidean system pQ. In the strong coupling phase (/? < (3 C ), monopoles form global and 
complicated clusters and make the vacuum complex bringing about repulsive forces among 
Dirac eigenvalues [U IU H]. This repulsive force among eigenvalues, which is observed as 
the Wigner distribution in the neighboring level spacings, forms the non-vanishing spectral 
density at the spectral origin, which is equivalent to the non-vanishing chiral condensate 
(iptp). On the other hand, in the weak coupling phase {fl > (3 C ), large monopole clusters 
vanish and the vacuum structure is much simpler, which leads to weaker repulsive forces 
among eigenvalues [TJ HJ [5] and results in the Poisson statistics in the level spacings. This 
weaker repulsive force is not so strong that Dirac eigenvalues can form non-zero spectral 
density at the origin, and the chiral condensate vanishes in this phase. 

We in Ref. PQ varied the couplings around the critical coupling f3 c employing isotropic 4D 
systems with the total volumes fixed, and investigated the natures of low-lying Dirac modes. 
In this case, the chirally symmetric vacuum realized at j3 > C is rather simple and pertur- 
bative. It is however known that the "finite temperature" transition still exhibits several 
nonperturbative features. For example, in the case of QCD, the formation of the strongly- 
coupled quark-gluon plasma phase (sQGP) [H [7J [5] just above the transition temperature 
was recently suggested and attracting many interests. Nonperturbative aspects above the 
transition temperature T c can be also observed in the Wilson loops. Whereas temporal Wil- 
son loops exhibit deconfinement feature above T c , spatial Wilson loops still show an area law 
(i.e. temporally deconfined and spatially confined), which hints that the vacuum structure 
is not simple but still complicated even above T c . The (3+l)D compact QED has similar 
properties, and the remaining nonperturbative structures can be also found in it. We inves- 
tigate the properties of Dirac eigenmodes, eigenvalues, and monopoles in quenched (3+l)D 
compact QED, and see what happens. 
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The organization of this paper is as follows: In Sec. [H we briefly show our formalism. 
Several properties of low-lying Dirac modes and eigenvalues are clarified in Sec. [31 Sec. Q] is 
devoted to the discussion based on numerical results. We summarize the paper in Sec. [5] 



2 Formalism 

We adopt the Wilson gauge action at /3=1.01 for gauge fields, 

/^^(l-cos^Or)), (1) 

and employ the overlap-Dirac operator [5J 110) , which is constructed as 

Hw 



D = p[l+ l5S gn(H w )}=p l+ l5 —^= , (2) 

VHw _ 

and realizes the exact chiral symmetry on a lattice [TT1 [T2] . Here, Hw = 75 (Dw — p) is the 
hermitian Wilson-Dirac operator defined with the standard Wilson-Dirac operator Dw- The 
"negative mass" p is chosen in the range of < p < 2, which we set 1.6 throughout this paper. 
We impose the periodic boundary conditions in all the spatial direction for the fermion fields, 
and the anti-periodic boundary condition is imposed on the temporal boundary. We compute 
lowest 50 eigenpairs at each "temperature" , with 48 gauge configurations. All the eigenvalues 
Ai a t of D, which lie on a circle with the radius of p in a complex plain, are stereographically 
projected onto the imaginary axis via Mobius transformation [13j . 

A= 1-AW2P (3) 

The spatial volumes are all fixed to 12 3 and we adopt 5 different temporal lengths, N t = 1/T 
= 4, 6, 8, 10, and 12. The "finite temperature" phase transition in compact QED was 
extensively investigated [T3] and the "transition temperature" at (3 = 1.01 was found to lie 
around N t ~ 6. (See also Ref. [IS].) We expect that our setup can cover the systems below 
and above the transition temperature, though finite volume effects might not be negligible 
and the phase transition temperature would be slightly modified. 



3 Low-lying Dirac eigenvalues and corresponding modes 

In Ref. PP, we found anti-correlations between (near-zero) Dirac modes and monopoles, 
which indicates that near-zero modes are "scattered" by monopoles, and it was conjectured 
that monopoles as impurities for near-zero modes are responsible for the "complexity" of 
a vacuum and control the level dynamics of low- lying Dirac eigenvalues pQ. Abundance of 
monopoles implies a complex vacuum and the level spacing distributions of Dirac eigenvalues 
obey the Wigner distribution. On the other hand, absence of monopoles makes a vacuum 
simple and the Poisson-like distribution appears. The Wigner (Poisson-like) distribution 
is the consequence of strong (weak) repulsive forces among Dirac eigenvalues, and strong 
(weak) repulsive forces are responsible for a non- vanishing (vanishing) spectral density at the 
spectral origin. The spectral density at the origin is finally related to chiral condensate {'tptp) 
via the Banks-Casher relation [3J. Chiral condensate (tfnp) simply reflects the complexity of a 
vacuum, which is brought about by monopoles' degrees of freedom. These results in Ref. pQ 
were all obtained in isotropic systems where the temporal length N t is equal to the spatial 
extent N s (N t — N s 3> l/T c ). Then, is this scenario the case also at "finite temperature" 
(Nt < N s )7 The answer would be no. Though the chirally symmetric vacuum in the weak 
coupling regime is actually simple pQ , we expect that the symmetric vacuum realized at high 
temperature is still complex. In fact, monopole density is not much reduced even at high 
temperature (small Nt) region, which can be found in Fig. [TJ 

We go further with the clarification of the chiral phase transition in compact QED. The 
keywords are (1) spectral densities at the spectral origin, (2) level spacing distributions of 
Dirac eigenvalues, (3) correlations between low-lying Dirac modes and monopoles, and (4) 
spatial distributions of low-lying Dirac modes. 
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Figure 1: The normalized monopole densities (total monopole lengths divided by total volumes) 
are plotted. The solid line is drawn at Nt ~ 7 for reference. 
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Table 1: The numbers of exact zero- modes found in 48 gauge configurations at each inverse 
temperature Nt are listed. The i-th row gives the number of configurations with 0, 1, 2 zero- 
mode (s), respectively. 



3.1 Spectral densities and level spacing distributions 

We show in Fig. [2] the histograms of Dirac eigenvalues at N t =12, 8, and 4. The numbers of 
N,=12 N t =08 N t =04 




Figure 2: The histograms H ev (X) of the eigenvalues A of the overlap-Dirac operator D are 
plotted for each A^(12,8,4). The horizontal axis denotes —iX. All the eigenvalues Ai a t lying on 
a circle in a complex plain are stereographically projected onto the imaginary axis via Mobius 
transformation, A = — ^ lat /o . 

zero modes found in 48 gauge configurations at each inverse temperature l/T = N t (4, 6, 8, 
10, 12) are listed in Table[T] The histograms around the spectral origin at N t =8, 12 are found 
to be rather flattened and have non-zero heights, while that at 7V t =4 exhibits the vanishing 
density. These densities are related to the non-vanishing (vanishing) chiral condensate at 
Nt > 6 ~ 1/T C (N t < 6 ~ 1/T" C ) via the Banks-Casher relation. It is remarkable that 
the chiral phase transition is almost accompanied by the deconfinement transition, which 
would occur around N t ~ 6. The vanishing density seen in the 12 3 x 4 system could be 
the consequence of some finite (spatial) volume effects. In order to ensure that it is not the 
case, we show the normalized histograms H CV (X) obtained with 12 3 x 4, 14 3 x 4, and 16 3 x 4 
lattices in Fig. [3Jleft]. All three histograms exhibit the same behavior, and we can neglect 
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Figure 3: Left: The histograms H CV (X) of the eigenvalues A of the overlap-Dirac operator 
D, which are obtained with 12 3 x 4, 14 3 x 4, and 16 3 x 4 lattices. Middle and Right: 
The distributions of lowest nonzero eigenvalues at N t = 12 and 4 obtained with v = gauge 
configurations. 

the finite volume effects in 12 3 x 4 system. The lowest-nonzero-eigenvalue distributions at 
Nt = 12 and 4 in v = sector are respectively plotted in Fig. [3(middle] and Fig. [3jright] , for 
the purpose of reference. 

As actually demonstrated in Refs. [JJ [TB] , the spectral density at the spectral origin is 
enhanced by the "repulsive force" among eigenvalues, and this force can be clarified with 
the neighboring level spacing distributions of Dirac eigenvalues. We show in Fig. |4] the 
unfolded neighboring level spacing distributions P\ a t(s) obtained at N t =12, 8, and 4. The 
solid lines denote the Wigner distribution function Pwig(s) = f|s 2 ex p(~^s 2 ), which is a 
good approximation of the original distribution obtained by the random matrix theory with 
chiral unitary ensemble, and the dashed lines the Poisson distribution function Pp i(s) = 
exp(— s). We again find the Wigner-Poisson transition, which was also found in Ref. [1] and 
is consistent with the spectral densities at the origin. 
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Figure 4: The unfolded nearest-neighbor level spacing distributions -Pi a t( s ) at Nt=12, 8, and 
4. The solid lines denote the Wigner distribution function fWig(s) = p-s 2 exp(— ^s 2 ) and the 
dashed lines the Poisson distribution function Pp i(s) = exp(— s). 



3.2 low- lying Dirac modes and monopoles 

We next investigate the correlations between low-lying modes and monopoles following 
Ref. [TJ. For this aim, we define and investigate the histogram ratios R^(p^) [JJ defined 
as 
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= ^iTTTT-- ( 4 ) 
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The eigenmode density p^ = \ip\(x)\ 2 here is the absolute squares of an eigenfunction 
^a(^) locally summed up over the spinor index. H^(p^) (HT on (p^)) denotes the histogram 
of p^ evaluated at all the sites (only on monopoles). This quantity R^,(p^) equals to 1, if 
there is no correlation between the spatial fluctuations of Dirac modes and monopoles. In the 
case when a positive (negative) correlation exists between the spatial fluctuations of Dirac 
modes and monopoles, R^(p^) > 1 at smaller (larger) p^ and i?^(p^,) < 1 at large (smaller) 
Pi/, hold. We show in Fig. [5] the histogram ratios for near-zero modes obtained at 2Vf=12, 8, 
and 4. Display ranges are chosen as < p^ < 0.0002 x j^- with intrinsic densities in mind. 
The universal anti-correlations are again observed. (We here omit the results on exact-zero 
modes |17[ 118] , since they show no apparent correlation with monopoles. This tendency was 
also found in our previous paper pQ. The present strategy is rather simple and we would 
need more sophisticated analyses to cast light on the properties of exact-zero modes. In the 
following sections, attentions will be paid mainly to near-zero Dirac modes.) 
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Figure 5: The histogram ratios R^j(p^) for near-zero modes at A^=12, 8, and 4 are plotted in 
the range of < p^ < 0.0002 x j^-. We draw a line at R 1 p(p 1 p) = 1 for reference. 

3.3 Spatial distributions of low-lying Dirac modes 

To clarify the spatial distributions of Dirac eigenmodes, it is convenient to extract the in- 
verse participation ratio (IPR) of each eigenmode. The eigenmode's distribution is partially 
encoded in the corresponding IPR: The IPR is unity when tp\{x) maximally spreads over 
the system, and equals to V in the case when ip\( x ) lives only on a single site, reflecting the 
spatial distribution of the eigenfunction ip\(x). 
The IPR /(A) of ipx(x) is defined as 

I(\) = vY / Pi?a(x) 2 , p IPR (x)^Y,\Mx)\ 2 - (5) 

x a 

Here, V denotes the system volume and ip\ (x) is the eigenfunction associated with an eigen- 
value A normalized as J2 X IV'aMI 2 = 1< The Greek alphabet a is the index for a spinor. The 
density /?ipr(x) is obtained by locally summing up the absolute square of each component 
of an eigenfunction ip\(x) only over its spinor index. 

We show the scatter plots of the IPRs of low-lying modes in Fig. [6] The IPRs of near-zero 
modes approximately range from 1 to 2 at N t =8 and 12. On the other hand, at N t =A, many 
of the IPRs of near-zero modes take much larger values indicating the strong localizations of 
near-zero modes above T c . Near-zero modes are localized and the overlaps among them get 
much smaller, which would be the source of the Poisson-like distribution in the neighboring 
level spacings of Dirac eigenvalues: Level repulsions would not occur without overlaps among 
eigenfunctions. Weaker repulsive forces among eigenvalues characterized by such a Poisson- 
like statistics lead to the vanishing eigenvalue density at the spectral origin, and the chiral 
symmetry is restored. The chiral phase transition across T c is now found to be essentially 
different from that across (3 C pQ. The Poisson-like statistics above (3 C , where near-zero Dirac 
modes are completely delocalized and are plane-wave like, is caused by the "simpleness" of 
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the vacuum. That above T c is however driven by the localization of near-zero modes. We 
expect nontrivial vacuum structures above T c . 
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Figure 6: The scatter plots of inverse participation ratios of low-lying Dirac modes at A^=12, 
8, and 4. The horizontal axis denotes —iX, the associated eigenvalue to each Dirac mode. 

We define another type of IPR Js(A;/u), which is nothing but the IPR evaluated in a 
3-dim subspace at x H = const., for further clarification of Dirac modes' distributions above 
T c . The subspace-IPR Is(X; a) is defined as, 

Is{X;^)=V s Y, PlFA&tf, pW^/^^EI^^I 2 / E \Mx)\ 2 - (6) 

x^— const, a a,x,x^— const. 

Here, V s denotes the volume of the 3-dim subspace. We evaluate the ratios, Is(X; i)/I(X)(i = 
1,2,3) and Is(X; 4)/J(A). In case an eigenfunction is simply extended along the temporal 
direction and spatially localized, 7s(A;4)//(A) is just unity and Is(X;i)/I(X)(i — 1,2,3) is 
expected to be less than 1, because in this case the eigenfunction has no structure along the 
time-direction and only the spatial structure is responsible for the total IPR. We show in 
Fig. [7] Is (A; 4) //(A) and Is(X;i)/I(X)(i = 1,2,3) obtained in 12 3 x 4 lattice. One can find 
that 7s(A;4)//(A) ~ 1 and Is(X;i)/I(X)(i — 1,2, 3) 1 hold, which indicates the vanishing 
temporal structure of low-lying Dirac modes at N t — 4. All the ratios, 7s(A;4)/J(A) and 
Is(X;i)/I(X)(i = 1,2,3), evaluated at 4 < N t < 12 are shown in Fig. [U The temporal 
structures of low- lying modes quickly vanish around Nt ~ 6 ~ 1/Tc. 
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Figure 7: The IPR-ratios, Is(\;i)/I(X)(i = 1,2,3) [left] and Js(A;4)/J(A) [right] obtained with 
the 12 3 x 4 lattice are plotted. 



4 Discussions 

Now that we have found that, at "high temperature" , low- lying (near-zero) Dirac modes lose 
their temporal structures and are spatially localized. We have two questions: (1) Why do the 
temporal structures of near- zero modes vanish? (2) Why are near-zero modes spatially (and 
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Figure 8: The IPR-ratios, Is(\;i)/I(\)(i = 1,2,3) [left] and Js(A;4)/I(A) [right] obtained with 
12 3 x 4 ~ 12 3 x 12 lattice are plotted. 



strongly) localized? As we have just found or as can be found in Ref. pQ, near-zero Dirac 
modes have anti-correlations with monopoles and then monopoles are surely responsible for 
the distributions of near-zero modes. These features may come from the metamorphosis of 
monopole world lines. 

Monopole world lines will take "static" configurations at high temperature; most of the 
world lines are twisted around the 4-dim torus along the temporal direction as is illustrated 
in Fig. [9jleft]. In such a case, near-zero modes can be freely extended along the temporal 
direction, since they do not encounter monopoles as impurities for near-zero modes, which is 
schematically illustrated in Fig. [9jmiddle] . This "clearing up" along the temporal direction 
may be the main cause for the vanishing temporal structures of near-zero modes. It could 
be also the origin of the stronger localization properties of near-zero modes: The clearing up 
leads to the irrelevance of the temporal direction and implies that the system dimensionality 
for near-zero Dirac modes is reduced approximately to three. The near-zero modes are then 
expected to get strongly localized due to this dimensional reduction. (Dimensionality also 
plays an essential role in the localization of wavefunctions. It is a celebrated fact that, in a 
system whose dimensionality is lower than two, any random potentials readily lead to the 
exponential localization of wavefunctions.) 

One may wonder which is the main origin of the effective dimensional reduction, clearing 
up along the temporal direction or shorter temporal extent N t ("thermal" effect). To answer 
this question, we perform exactly the same analyses in an artificially constructed isotropic 
12 3 x 12 system. We define link variables C/^ 2 '(x, t) for this new system with {/^(x, t) in the 
original 12 3 x 4 lattice. U^\k-, t) is defined as 

Ujp(x,t) = U^t mod 4) (0<*<11). (7) 

In other words, we construct an isotropic 12 3 x 12 system by piling up three 12 3 x 4 systems. 
(See Fig. [Qjright].) Due to the periodicity of U(l) gauge fields, monopole configurations 
are the same in these two systems. Though this newly constructed lattice is an isotropic 
system, monopole world lines exhibit static configurations. We show in Fig.QlDlthe histogram 
ratio R^(p^j) for near-zero modes, the scatter plot of inverse participation ratios of low- lying 
Dirac modes, the unfolded nearest-neighbor level spacing distributions -Piat(s), the IPR-ratio 
Is(\;i)/I(X)(i = 1,2,3), and the IPR-ratio Is(X; 4)/7(A). Surprisingly, all the quantities 
essentially remain the same, which indicates that the properties of low-lying modes do not 
depend on the temporal lattice extent at all, and that the temporal dimension is actually 
irrelevant. 

The chiral phase transition in compact QED system is controlled by the metamorphosis 
of monopole world lines: When a system is covered by large monopole clusters, the system is 
complex p] and the spectra of Dirac operators exhibit the Wigner distribution, which implies 
strong repulsions among Dirac eigenvalues and gives rise to the non- vanishing spectral density 
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Figure 9: Left: Monopoles take rather static configurations with small Nt (at high temperature). 
Middle: Near-zero Dirac modes are localized avoiding monopoles. In the presence of static 
monopoles, near-zero modes would be spatially localized and temporally delocalized. Right: 
We artificially construct 12 3 x 12 system duplicating 12 3 x 4 system. Due to the periodicity of 
gauge fields, monopole configurations on this artificial 12 3 x 12 system are the same as those in 
the original 12 3 x 4 system. 

at the origin ((tptp) ^ 0). On the other hand, at high temperature, monopole world lines take 
static configurations and the system undergoes the clearing up along the temporal direction, 
which would lead to the effective dimensional reduction for near-zero modes and to the 
exponential localization of the modes. The exponential localization of Dirac modes causes 
weaker repulsive forces among eigenvalues exhibiting the Poisson-like distribution of Dirac 
spectra. Such weak repulsive forces cannot form non-zero spectral density at the spectral 
origin ({ipip) = 0). The metamorphosis of monopole world lines is also responsible for the 
deconfinement phase transition. In a system where large monopole clusters cover the entire 
volume, it is in confinement phase. When monopole world lines take static configurations, 
the system breaks out from the confinement phase to the deconfinement phase. Such static 
monopoles still lead to the spatial confinement, which is the remnant of the confinement and 
gives rise to remaining nonperturbative phenomena. 

These two phase transitions are then considered to be both controlled by the monopole 
dynamics. The coincidence of two "transition temperatures" could be naturally understood 
from this viewpoint. At least, if the metamorphosis of monopole world lines in the vicinity 
of the phase transition point is quick enough, the transition temperatures would be similar 
in magnitude. The clarification of Polyakov loops in terms of Dirac spectraQlJ] could also 
give us deeper understanding of the possible relationship between confinement and chiral 
symmetry breaking. 

5 Summary 

We have studied the properties of low-lying Dirac modes in quenched compact QED at 
£=1.01, employing 12 3 x N t lattices with N t = l/T = 4, 6, 8, 10, 12. The overlap formalism 
has been adopted for the fermion action. 

We have found several features worth noting: 

• The nearest-neighbor level spacing distribution of Dirac eigenvalues coincides with the 
Wigner distribution at Nt = 12 (confinement phase), and it gradually changes to the 
Poisson distribution as Nt is decreased, which is consistent with the spectral density 
at the spectral origin. 

• Near-zero modes exhibits stronger localization at high temperature, which is the origin 
of the Poisson-like distribution of Dirac spectra. 

• Near-zero modes lose their temporal structures above T c . 

• Near-zero modes have been found to have universal anti-correlations again with monopole 
world lines below and above the critical temperature. 

The chiral phase transition at "high temperature" in compact QED seems to be controlled by 
the metamorphosis of monopole world lines. Monopole world lines take static configurations 
at high-T system and the system undergoes the clearing up along the temporal direction, 
which would lead to the effective dimensional reduction for near-zero modes. This reduction 
is responsible for the exponential localization of the modes, which causes weaker repulsive 
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Figure 10: Several quantities obtained with the artificially constructed 12 3 x 12 lattice are 
plotted; upper left: The histogram ratio R^(p^) for near-zero modes, upper middle: the 
scatter plot of inverse participation ratios of low- lying Dirac modes, upper right: the unfolded 
nearest-neighbor level spacing distributions -Piat(s), lower left: the IPR-ratio Is(\;i)/I(X)(i = 
1,2,3), lower right: the IPR-ratio Is(X; 4)/T(A). Surprisingly, all these results remain almost 
the same as those in 12 3 x 4 system. 

forces among eigenvalues and gives rise to the Poisson-like distribution of Dirac spectra. Due 
to such weak repulsive forces, the spectral density at the spectral origin vanishes ({ipip) = 0). 
Taking into account that the metamorphosis of monopole world lines is also responsible for 
the deconfincment transition, we conjecture from the microscopic viewpoint that the chiral 
and confinement/deconfinement transitions in compact QED are both induced by a single 
origin, monopole's dynamics. 
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